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Abstract. Finite-dimensional spaces which are invariant under the action of the 
Hamiltonian of the BCn Inozemtsev model are introduced, and it is shown that 
commuting operators of conserved quantities also preserve the finite-dimensional 
spaces. 

The relationship between the finite-dimensional spaces of the BCn Inozemtsev 
models and the theta-type invariant spaces of the BCn Ruijsenaars models is clar- 
ified. The degeneration of the BCn Inozemtsev models and the correspondence of 
their invariant spaces are considered. 



1. Introduction 

In 0, Inozemtsev proposed a iV-particle quantum mechanics model, which is 
called the BCn Inozemtsev model. This is a generalization of the Calogero-Moser- 



Sutherland model or the Olshanetsky-Perelomov model . 

The BCn Inozemtsev model is quantum completely integrable. Here, the quantum 
complete integrability means that there exists N algebraically independent mutually 
commuting operators (higher commuting Hamiltonians) which commute with the 
Hamiltonian of the model. It is a quantum version of Liouville's integrability. For the 



BCn Inozemtsev model, Oshima |13] described the commuting operators explicitly. 
Note that the BCn Inozemtsev model is a universal completely integrable model of 
quantum mechanics with Bn symmetry, which follows from the classification due to 
Ochiai, Oshima and Sekiguchi (JTTj, |TJJ. 

On the other hand, Finkel, Gomez-Ullate, Gonzalez-Lopez, Rodriguez and Zh- 
danov studied quasi-exactly solvable models in 0, §]. They found several quasi- 
exactly solvable many-body systems. Although they did not use the phrase 11 BCn 
Inozemtsev model", they essentially found that the BCn Inozemtsev model is quasi- 
exactly solvable, i.e. the Hamiltonian of the BCn Inozemtsev model preserves some 
finite-dimensional space which is spanned by some symmetric "monomials" . 

In this paper, we link the quasi-exact solvability with the quantum complete inte- 
grability. More precisely, we show that the commuting operators (higher Hamiltoni- 
ans) of the BCn Inozemtsev model also preserve the finite-dimensional space, which 
has been appeared in the context of quasi-exact solvability. 

On the finite-dimensional space, joint eigenvalues and eigenfunctions of the com- 
muting operators are determined by algebraic calculations. In this sense, the model 
would be solved partially. Note that the phrase "quasi-exact solvability" is used in 
these situations (see pl|j ). 



The spectral problem of quantum mechanics is generally considered in a Hilbert 
space, and the Hilbert space is often taken as a square- integrable space [L 2 space). 
Therefore it would be important to consider the relationship between the Hilbert 
space {I? space) of the BCn Inozemtsev model and the finite-dimensional space which 
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appears in the context of the quasi-exact solvability. In this paper, we determine the 
condition for that the finite-dimensional space lies in the L 2 space. 

There are other models which are concerned with the quasi-exactly solvability. 
Ruijsenaars-type models are introduced for arbitrary root systems including the BCn 
cases, which are difference analogues of Inozemtsev (or Calogero-Moser-Sutherland) 
models (see JO)], |T|, Hikami and Komori |^, 0, |8|, [J finally constructed higher com- 
muting operators using algebra, and found an invariant subspace spanned by theta 



functions. In ||16|| , Sasaki and Takasaki considered degenerate Inozemtsev models and 
their quasi-exact solvability. 

In this paper, correspondence between the BCn Ruijsenaars model and the BCn 
Inozemtsev model is considered. We will observe that the invariant subspace spanned 
by theta functions for the BCn Ruijsenaars model corresponds to the invariant space 
related to the quasi-exact solvability for the BCn Inozemtsev model. 

To obtain the degenerate BCn Inozemtsev model from the (elliptic) BCn Inozemt- 
sev model, a certain trigonometric limit is considered. It is shown that the finite- 
dimensional invariant spaces for the BCn Inozemtsev model tend to the invariant 
spaces of degenerate BCn Inozemtsev model which were introduced by Sasaki and 
Takasaki. It would be important to consider degeneration of models, because it would 
be helpful in understanding several integrable models and the relationship among 
them. 

This paper is organized as follows. In section [| the BCn Inozemtsev model and its 
finite-dimensional invariant spaces are introduced. In section |3.1], higher commuting 



operators of the BCn Inozemtsev model are introduced. In section p.2|, it is shown 



that higher commuting operators also preserve the finite-dimensional invariant spaces. 
In section |3.3| , the relationship between the finite-dimensional invariant spaces and 
the L 2 space is considered. If the coupling constants are integers, the model may 
have some special features. In section |^|, we consider this case. In section f|, corre- 
spondence between the theta-type invariant spaces for the BCn Ruijsenaars model 
and the invariant spaces which is related to the quasi-exact solvability for the BCn 
Inozemtsev model are investigated. In section [|, we consider the degeneration of the 
BCn Inozemtsev models and the limit of invariant spaces. 
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Oshima, Ryu Sasaki, and Kanehisa Takasaki for fruitful discussions. He is partially 
supported by the Grant-in-Aid for Scientific Research (No. 13740021) from the Japan 
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2. BC N Inozemtsev model and its invariant space 

The BCn Inozemtsev model is a system of quantum mechanics with iV-particles 
whose Hamiltonian is given by 

N d 2 

(2.1) # = _^_ + 2/(Z + l) {p(xj - x k ) + p{ Xj + x k )) 

j=l 3 l<j<k<N 
N 3 

+ Yl Yl + + w <)> 

j=l i=0 
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where p(x) is the Weierstrass' p-function with periods (l,r) (see ( |A.1| )), uj = 0, 
lji = |, u 2 = — 2 -j^-, = | are half periods, and Z and l, t (i — 0, 1, 2, 3) are coupling 
constants. 




Let a, 6j (z = 0, 1, 2, 3) be numbers which satisfy a G {— Z, Z + l} and bi G {— 
(z = 0, 1, 2, 3). Set ^ = p(xj-) (1 < j < N) and 

N 3 

(2.2) *(*)= n (%--^) a nn^- e ^> 

l<j<k<N j=l i=l 

H = $(z)- 1 oH o$(z), 

where = p(cu;), (z = 1, 2, 3). 

By applying formulae ( |A.3j ), it is shown directly that the operator H admits the 
following expression: 

(2.3) 

/ N ( d 2 

H = - \ J2^( z j - ~ e 2 )(*j - e 3 ) 

2a 26i + i 26 2 + | 

+ > + " + " 

z i ~ z * z i - e i 

- 4 UN - l)a - b + 61 + 6 2 + h + ^\ ((N - I) a + b Q + h + b 2 + b 3 ) 

+ 4jV((&! + b 2 ) 2 e 3 + (b x + b 3 ) 2 e 2 + (6a + 63)^1) 

- 4iV(N - l)a(e x &i + e 2 6 2 + e 3 6 3 ). 

Proposition 2.1. (%) LeZ P s y m be the space of symmetric polynomials in variables 
z u z 2 ,...,z N . Then H ■ P s ^ a C P sym . 

(ii) Let a, bi (i = 0, 1, 2, 3) be numbers which satisfy a G {— Z, Z+l} and 6j G { — §, ^p} 
(i = 0, 1, 2, 3). Assume that d = —((N — l)a + 6 + &i + b 2 + b 3 ) is a non-negative 
integer. Let Vd be the vector space spanned by monomials z^z™ 2 . . .z™ N such that 
mi G {0, 1, ... , d} for all i, and V° ym = V d f] P^ m . Then H ■ V* ym C V* ym . 

Proof. Let / G P sym . From Q2.3| ), the function Hf is a symmetric rational function 
and it may have poles only along Zj — Zk — of degree at most 1. If Hf has a pole 
along Zj — Zk = (j 7^ fc) of degree 1, it contradicts to the symmetry of Hf on the 
variables Zj and z^. Hence the function Hf does not have poles and we obtain (i). 

Let Ai, . . . , Xn be non-negative integers such that Ai > \ 2 > ■ ■ ■ > \n > 0. 

Set P{z) = 4(z — ei)(z — e 2 )(z — e 3 ) and L = Ai. Let Z and I' be integers such that 
< Z, V < L. Then 

(2-4) ( 7~~-~T~^7~ + T^T^l (44' + 4'4) 




(zj-z k )dzj (z k -Zj)dz k/ 

ALS^izf+'zi + z i ; z L k +l ) + AL5 v>L {z l jZ ^ x + + (#1;, 



where is the Kronecker's delta and the term (#1) is a linear combination of 
monomials z^z^ such that < t,t' < L. 
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Let A = (Ai, A 2 , • • • , Xn) and A + = (Ai + 1, A 2 , . . . , Xn). By summing up the equality 
on j and k (1 < j < k < N), we obtain 



N P( \ P) 

where m M = T l { mi ,..., mN )e& N -n z T lz T 2 ■■■z% N for /i = (/ii, . . . ,^ N ), S N is the sym- 
metric group of N letters and the term (#2) is a linear combination of symmetric 
monomials m^...^) such that < /xjv < • • • < A*i < 
Hence we obtain 

(2.6) #m A = -4(L - d - 26 + - d)m x+ + (# 3 ), 

where the term (#3) is a linear combination of symmetric monomials rn.^...^) such 
that < hn < • • • < A*i < £■ 

If L < d then all elements in the term (#3) lie in V^ sym . If L < d then m A + G V^ ym , 
and if L = d then the coefficient of m\+ on the right-hand side of ( [2.6|) vanishes. 
Hence if L < <i then ifm A G V^ sym . 

Therefore we obtain (ii). □ 



Remark Proposition |2.1| was essentially obtained in P) by a different method. □ 



3. Commuting operators and invariant subspaces 

3.1. Commuting operators. 

It is known that the BC^ Inozemtsev model is completely integrable, i.e. there 
exists N algebraically independent mutually commuting operators (higher commuting 



Hamiltonians) which commute with the Hamiltonian. In |13| , Oshima gave explicit 
forms of the commuting operators. Now we pick up some results obtained by Oshima. 

Let W(Bn) be the Weyl group of type Bn, i.e. the group of the coordinate trans- 
formations 



(3.1) (xi, . . .,x N ) h-> (eiX CT( i), . . . , e N x a ( N )) 



of Mr, where a G &n (the symmetric group) and e\ = ±1, . . . , e N = ±1. Let W(D 
be a subgroup of W(B^) which consist of transformations ( |3.ip with a condition 

Ylf=i *j = 1. For w G W(B N ) we define e(w) = { \ ™ | JJ^j . 
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Let us consider the operators which commute with the Hamiltonian of the BC^ 
Inozemtsev model fl2TT|). For i — 0, 1, 2, 3, we set 

(3.2) S {1} ... :k} = ^ w(p(x 1 -x 2 )p(x 2 -x 3 )...p(x k ^ 1 -x k )), 

weW(B k ) 

S tl-,k} = E w(p(xi +Ui)p(xi -x 2 )p(x 2 - x 3 ) . . . p(x k -i -x k )), 

w&W{B k ) 

T° { i,..,k}= E Hr 1 ^-!)^...^, 

/ 1 U-U/ M ={l,...,fe} 
7iU-U7 M ={l,...,fc} 

3 

T{i,...M = + l))^ 1 E ^^ { °£,*P 

where the sum ^^n.-n/ ={i runs over a ^ different partitions of {1, . . . , A;}. For 
example 

TO 1 rpo _ Q rpo _ Q GO 

_ ' {1} ~~ '"'i 1 }' {1.2} ~~ ^{1.2} ~~ >J{l}>->{2}, 

^{1,2,3} = ^{1,2,3} - 5'{1}'S'{2,3} - S{ 2 }S{ li3 } - 5*13} S{ 1>2 } + 25{i}S'{2}5'{3 } . 

Set 

(3.3) A {li ... ifc} = 2*Vk-2')\ S e ( w ) w (p( x i -x 2 )p(x 3 -x 4 ) ... 

o<i<[|] ^ wew(B k ) 

. . d d d 

■<P( X 2j-\ ~ X 2 j)- 



dx 2j+1 dx 2j+2 dx k j 

Q{i,...,k}= E T h--- T i», 
/in..-u/ M 6{i,...,fc} 

T w ({i,... t k}) = w(T {1 ^ k} ), A w{{lt _ jk}) = iw(A{i ) .„ )fe }), for u> G 6jv, 
where [|] represents the maximum integer not greater than |. 



Proposition 3.1. (JT3[ theorem 7.2, remark 7.4)- Set 

N N 

(3 - 4) p ^ = EE,-u,_,-w/v_,h E E 



z!(j -i)!(JV- j)! 
i=fc j=« w ' ^ •> ' we& N hu-ui k ={i,...,i} 



w 



{(-1(1 + l)r k 2- k Tl . . .T? k q {i+1 ^ y Al j+1 _ N} ) 
for k = 0, . . . , N — 1. Then the operators Pj (1 < j < N) are W (B n) -invariant and 
(3.5) [Pj,P k ] = 



for 1 < j, k < N. The Hamiltonian H ( 2.J\) is a linear combination of Pi and 1, i.e. 
H = AP\ + B for some constants A, B. 
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3.2. Commuting operators and invariant subspaces. 

We change variables by Zj — p(xj) (j — 1, . . . , N) and set 

(3.6) P k = <S>{z)- l oP k a<S>{z), (k = l,...,N) 
where Q(z) is defined in ( |2.2| ). From ( |3.5| ), we obtain 

(3.7) [Pj,h] = 
for 1 < j, k < N. 

Proposition 3.2. The operators P k (k — 1, . . . , N) admit the following expansion: 



(3.8) P k = 



0<ii,...,i N <2 
l"*JV— 2fc 



dzij \dz N 



Here, the operators Pk are symmetric in Zi, . . . , zjy and the coefficients Ci lt ,,, t i N (z) are 
rational functions in z\, . . . ,z^ which may have poles only along Zj — = and 
z h ~ z j2 = (1 < i < 3, 1 < j,ji,j 2 < N, ji ^ j 2 ). 



Proof. From (|3.4j) , the operator Pk admits the expansion 

Q \ il / Q \iN 



(3.9) P k = d ix,-,iN( x h---,X N ) [^-J 



0<i_ 

hijv<2fc 



such that di lr .. s i N (xi, . . . ,xn) is doubly periodic for each variable xi, . . . ,xn- From 
the iy(Piv)-symmetry of the operator Pk, we obtain di lt _ i i N (xi, . . . , —Xj, . . . ,xn) = 
(-l)^d iu ,^ iN (x 1 , . . . ,Xj, . . . ,x N ) for j = 1, . . . , N. 

Since p'(x) is an odd doubly periodic function, the function d\ iN [x\, . . . , x?q) 
defined by ^,...,^(^1, ...,x N )= d? u ... )ijv (zi, ...,x N ) Ylj. i:j=1 p'(xj) is even doubly pe- 
riodic in each Xj (j = 1, . . . , N). 

Now, change variables Zj = p(xj) (j = 1, . . . , N) and write 



d \ n ( d x %N 

(3.10) P k = V, dn^i N {x x , . . . ,X N ) 



0<i lt ...,i N <2 
ilH MiV— 2fc 



dz 1 J \dz N 



Since p'(x) 2 and p"(x) are even doubly periodic, the function d ilt ^^ iN (xi, . . . ,xjy) 
is even doubly periodic in each Xj (j = 1,...,N). It is known that even doubly 
periodic function in x is expressed as F(p(x)) with some rational function F(z). 
Hence if we set Zj = p(xj) (j = 1, . . . ,N), the functions di u _ t i N (xi, . . . ,Xjy) admit 
the expression di u _ t i N (xi, . . . , Xjv) = Cii,...,i N (p(%i), ■ ■ ■ , p(xn)) with some rational 
function Ci u ,„ t i N (zi, . . . , Zn)- From equality ( |3.4| ) and formulae ( |A.3| ), it is shown 
that the coefficients Ci lt ...,i N (zi, . . . ,z^) may have poles only along Zj — e, = and 
z h ~ z h = (1 < i < 3, 1 < jjijs < N, ji ^ j 2 ). 

Next we consider the coefficients of the operators Pk defined by ( |3.6| ). From prop- 
erties of the functions (k ly „ : i N (zi, . . . ,Zn), it is shown that the coefficients Ci lt ,„ t i N (z) 
in ( ^.8| ) are rational functions in variables z\, . . . , zn which may have poles only along 
Zj - ei = and z jl - z h = (1 < i < 3, 1 < j,ji,j 2 < N, ji ^ j 2 ). 

Since the operators Pk and the function $(2;) are symmetric in z\,...,zn, the 
operators Pk are symmetric in z\, . . . , z^. □ 
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Theorem 3.3. (i) Let P s y m is the space of symmetric polynomials in variables 
z x ,...,z N , then P k ■ P sym C P sym for k = 1, 2, . . . , N 

(ii) Let a, bi (i = 0, 1, 2, 3) be numbers which satisfy a G {— /, 1+1} and &j G { — 2^} 
(z = 0, 1, 2, 3). Assume that d = —((N — l)a + & + b\ + 62 + ^3) a non-negative inte- 
ger. Then Pk ■ V^ ym C V^ sym for k — 1, 2, . . . , N, where V^ ym is the finite- dimensional 
space defined in proposition W. 



Proof. The coefficients Ci lt ... t i N (z) in ( |3.8| ) are rational functions which may have poles 
only along zj - e { = 0, - ^ a = (1 < i < 3, 1 < J, Ji, J 2 < iV, ji ^ j 2 ). 

Let us fix ji and j'2 that satisfy j\ 7^ j2- Let p be the maximal number of degrees 
of a pole along z^ — Zj 2 = of functions (ki,...,i N {z) f° r possible all ii, - - - , *jv- If 
f(z) G P sym then the function P k f(z) has a pole at most degree p along 2^ — Zj 2 = 0. 
Let /(z) be an element of P sym such that the function P k f(z) has a pole of maximum 
degree along Zj 1 — Zj 2 = 0. We denote the degree by p' . Then the function P k Hf(z) 
has a pole of degree at most p' along Zj x — Zj 2 = 0, because Hf(z) G P sym . On the 
other hand, it is shown that the function HP k f(z) has a pole of degree p' + 2 along 
^ - 2j2 = if p' ^ 0, 1 - 2a. Since P k Hf(z) = HP k f{z), if 1 - 2a G" Z then p' 
must be equal to zero. Hence the function P k f(z) is holomorphic along — z j2 = 
if 1 — 2a (jL Z. By a continuity argument in a, it is shown that the function P k f(z) is 
holomorphic along Zj ± — Zj 2 = for all a. 

Similarly if f(z) G P sym then the function P k f(z) is holomorphic along Zj — = 

and - z j2 = (1 < i < 3, 1 < j, j 1; j 2 < ^ ^ j 2 ). Hence P fe /(z) G P sym . 
Next we prove (ii). 

From expression (|3.8| ), there exists p G Z> such that P k f(z) C V^ sy ™ for all 
f(z) G V^ sym . Let /(z) be an element of V^ sym such that the degree of P k f(z) is 
maximum. We denote the degree by d + p'. Then the function P k Hf(z) has a degree 
at most d + p', because Hf(z) G V^ sym . On the other hand it is shown that the 
function HP k f(z) has a degree d + p' + 1 if p 1 ^ 0, 26 — 1/2. From the commutativity 
of P and Pfc, if 26 — 1/2 ^ Z then p' must be equal to zero. Thus P k f(z) G V^ sym . 
By a continuity argument, we can remove the condition 2&o — 1/2 G" Z. 

Hence we obtain(ii). □ 

In summary, we established that the higher commuting Hamiltonians also preserve 
the space related to the quasi-exact solvablity in theorem [O] (ii). 

3.3. Relationship to the L 2 space. 

Assume IJ0J1 £ ^>o an d I2J3 £ K in this subsection. 

The invariant space V^ sym of the PCV Inozemtsev model is defined for each a G 
{-/, / + 1} and bi G {-§, ^} (z = 0, 1, 2, 3) with the condition d = -((N - l)a + 
bo + bi + b 2 + 63) G Z> . For these numbers a, 60, bi, 62, &3, define 

(3.11) W s / m = {$(p(x x ), . . . , p(x N ))f(p( Xl ), . . . , p(a*r))|/(*i, - - - , z N ) G V d sym }, 

where the function . . . , z N ) = ni<i<fc<iv( z i ~ z kY IlJLi IlLifo ~ e i) bt was 

defined in 



From relations (|2.2j ), ( p.6|) , and Theorem |3.3j , the following proposition is shown 
immediately: 
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Proposition 3.4. Assume that d = —((N — l)a + 60 + b\ + b 2 + ^3) is a non-negative 
integer. Then H ■ W* ym C W f m and P k ■ W* ym C W f m for k = 1, 2, . . . , N. 



in 

1 



We look into the condition that the space W^ ym lies in L 2 space. If f(x) G W^ y 
then ~ \xj—Xk\ a (resp. \f(x)\ ~ \xj+Xk\ a ) as \xj— Xk\ — > (resp. — > 0) 

(j 7^ fc)? 1/(^)1 ~ |xj| 2fe ° as — > 0, and ~ \xj — || 2bl as |rcj — \\ —> 0. Hence 

if a = / + 1, 60 = ^2^", and 61 = then the function /(x) G H^ ym is locally square- 
integrable on \xj — x^\, \xj + Xk\, \%j — || < e for sufficiently small e. Combining 
with the periodicity of f(x) G W^ ym , we obtain J" 0<a;i <-...<,,. <i \f(x)\ 2 dxi . . . dxN < 00. 
Hence the following proposition is shown: 

Proposition 3.5. Let b { G {-§, ^±1} (i = 2,3). 1/(2 = -({N - + 1) + ^ + 
1 + ^2 + ^3) £ Z> then every function in W^ ym is square-integrable on the domain 
< x\ < ■ ■ ■ < x N < 1. 

In the case d = -((N — l)(l + l) + l -^ x + l + b 2 + b 3 ) G Z> , some eigenvalues of the 
commuting Hamiltonians on the Hilbert space (L 2 space) appear as the eigenvalues on 
the subspace W^ ym . Hence some eigenvalues on the Hilbert space would be obtained 
explicitly, because the eigenvalues in the finite-dimensional space are got by algebraic 
calculations. Note that the case BC\ was done in [pl|, and Gomez-Ullate, Gonzalez- 
Lopez, and Rodriguez considered the relationship between the L 2 space and the space 
related to the quasi-exact solvablity for some special cases in [|J. 

As an aside, the joint eigenvalues of the trigonometric BC^ Calogero-Sutherland 
model are already known and their expression is simple. Distributions of eigenvalues 
in L 2 fl W^ ym will be detected by considering the trigonometric limit r — > \J — loo 
while fixing coupling constants l,l ,li,l 2 ,l 3 . 

3.4. The case I, l , l\, l 2 , h G Z> . 

Let us consider the case /, Iq, l\, l 2 , h G Z>o. In this case, the Hamiltonian and the 
higher commuting Hamiltonians preserve several finite-dimensional spaces of elliptic 
functions. 

The invariant space W^ ym is defined for each a G {—1,1 + 1} and foj G { — 2^} 
(i = 0, 1, 2, 3) with the condition d = -((N - l)a + b + bi + b 2 + 6 3 ) G Z> . 

For each I, l , l\, l 2 , l 3 G Z> , there are eight possible sets of (a, 60, 61, b 2 , 63) for each 
set the invariant space W^ ym is defined, if N > 2. If N — 1 then there are four 
possible sets of (b , bi, b 2 , b 3 ). For example, if iV > 2, I ^> l , / 1; l 2 , Z 3 and (N — 1)1 + 
lo + h + l 2 + I3 G 2Z >0 , then the cases 

{/_/ ^0 d± dil z±L\ (—] b±l h+1 ^2+1 h±l\ 
\ b i 2 ' 2 ' 2 ' 2 /' \ ' 2 ' 2 ' 2 ' 2 /' 
( — ] ^0 3\ h+1 h±±\ (_] 1o+l -h -h \ 

_ ' XiM XiAi) \Jx A ) 

\ t ' 2 ' 2 ' 2 ' 2 /' V ' 2 ' 2 ' 2 ' 2 /' 
(_j ^lo h±l h±l (1 h+l A ^ Z 3 +l \ 

V l i 2 ' 2 ' 2 ' 2 /' v ' 2 ' 2 ' 2 ' 2 /' 

are permitted. 

By a straightforward calculation, the dimension of direct sum of spaces W sym of 
elliptic functions can be calculated. For the case N — 1, the dimension is 



/ . 



2fc + 1, I is even and k + k 3 > „ . 

/ — 2k 3 + 1, I is even and k + k 3 < |; 

2A; + 1, f is odd and k > 

1 + 2, I is odd and k < 
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where / = l + l\ + l 2 + h, k = max(Z , ^1, h, h), an d k 3 = min(/ , h, h, h) ( see also 



T3, H3). For the case N = 2, the dimension is (21 + l) 2 + J2 i=0 Wi + l ) for a11 the 
case l,lo,h,h,h G Z> . Since the dimension is directly related to the genus of the 
spectral curve for the case N = 1 |I9 |, the dimension for the case N > 2 might also 
play important roles. 



4. RUIJSENAARS MODELS AND INOZEMTSEV MODELS 

In [[L5[], Ruijsenaars introduced a relativistic version of the Calogero-Moser-Sutherland 
model, which is called the Ruijsenaars model of type An these days. In van Diejen 
introduced the BCn Ruijsenaars-type model which has ten parameters (k; //, z/ , z/ , ^1, Pi, ^2, ^2, ^3, P3), 
and Hikami and Komori constructed higher commuting operators by use of root al- 
gebra in m [7|, H, that ensures the integrability. 

The lowest operator of the BCn (or -4^) Ruijsenaars model is given as follows: 



(4.1) 



N 



3=1 



N 

E 



n 

fc=i 

3 

n 

r=0 
( N 



\ 



gl (Xj -X k - fi) 61 (Xj +x k - jj) 

9x(xj-x h ) 6i(xj + x k ) ^ 

9 r+ i(Xj — Z/ r ) r+ i(Xj + k/2 — v r ) 

9 r +i{xj) e r+1 (xj + K/2) 

OxiXj + X k + /i) ^(Xj - Xfc + /i) 



n 



I 



n 



p=0 



9 r+1 (xj + z/ r ) 9 r+ i(xj — k/2 + z/ r ) 
1 6'. r+ i(x i ) 9 r+ i(xj - k/2) 

2 



tj(-K) 



71 



Jj0 P+1 («/2 + 



vr=0 



V 

n 

V7"=l 



9 p+1 (xj - k/2 - fi) 9 p+1 (xj + k/2 + /i ) 
d P +i{xj - k/2) 9 p+1 (xj + k/2) 



Here 9j(x) (j = 1,2,3,4) is the Jacobi theta function (see (|A.4|) ) and U(k) is a 
translation of the variable Xi by k, i.e. ti(n)f(xi, . . . , x i: . . . xn) = f(xi, ■ ■ ■ , %i + 
k, ...Xjsr). n r (r = 0,1,2,3) denotes the permutation tt = id, tti = (01)(23), 7r 2 = 

i, k = j. 

Hikami and Komori showed that if k = (2(N — l)/i + ^L (^i + ^)) I K e 2Z> 

then the operator Yi and higher commuting operators preserve the space of level k 

(2) 

theta functions of type A 2N . They proved it using root algebra. Their presentation of 
the invariant subspace would be technical for non-experts. In this section we describe 
them plainly. 
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The space of level k theta functions is defined as follows: 

{holomorphic, W (B ^-invariant 
f:C N ^C f(x + n) = f(x), ' (VneZ N ) 
f[x + nr) = /(a;)e-W=Tfc(W«)+N™W2) 

where (x\y) = J2?=i x iVi for x = { x i, • • • , x n) e C N and y = (y h . . . , y N ) G C N . A 
function f(xi, . . . , xn) is iy(-BAr)-invariant if and only if the relations /(x CT (i), . . . , x a ^)) 
f(xt, ...,x N ) for Va G 6at and /(a^, . . . , x h . . . , x N ) = f(x lt -x h ...,x N ) for 

Vi G {1, . . . , iV} are satisfied. For the case iV = 1, we obtain dimTh^^ 1 ^ = I + 1 for 
/ G Z> . Let flWfc), • • • , 6 {l+1) {x) be a basis of Th^ {Bl) . Then the space T/i^ (Bjv) is 
spanned by functions J2aee N ° {kl) ( x a(i)) ■ ■ ■ {kN) ( x a(N)) (1 < h < ■ ■ ■ < k N < I + 1). 
Therefore, dim TTi^" 8 ^ = f° r ^ e ^>o- 

Proposition 4.1. (c./. §) J/ Jfe = (2(iV - + ELo(^ + ^)) / K G 2Z >o i/ien t he 
operator Y\ preserve the space Tli^^ 3 ^. 

Proof. Let fix) G Th^^ BN \ Then the function Yif(x) is W^(-Bjv)-invariant. From 
the quasi-periodicity of 9i(x) (i = 0,1,2,3) (see ( |A.5| )) and f(x) (see (|4.2|)), the 
function Yif(x) has a quasi-periodicity as condition (|4.2|) when k = (2(N — l)fi + 



Si=o(^ + U i))/ K - Hence if we show the holomorphy of the function Y\f(x) on C N , 
we have Y\f{x) G TIi^^ Bn \ Thus it is sufficient to show that the residues of the 
function Y\f(x) at Xj — x^ = 0, Xj + x k = (1 < j ^ k < N), and Xj = 0, 1/2, (1 + 
r)/2,r/2,±«/2, 1/2±k/2, (1 + t)/2±«/2,t/2±k/2 (1 < j < AT) are zero. These are 



shown directly by using the quasi-periodicity of f(x) ( |4.2| ) and ^(x) (z = 0,1,2,3) 



( |A.5| ). Note that we rely on the condition k G 2Z in this step. □ 



Let us consider the non-relativistic (difference-differential) limit of the Ruijsenaars 
model. It is known that the Inozemtsev model appear by this limit. Now we will 
exhibit it explicitly. 

Let a = —/i/k, b = — (j/i + V\)/2k, b x = -{y 2 + v 2 )/2k, b 2 = —(u 3 + u s )/2k, 
h = -(^o + v )/2k and 
(4.3) 

JV 

Q(x) = H {e 1 {x j -x k )e 1 {x j + x k )) a l[e 1 {x j )^e 2 {x j )^e3{x j )^e (x j )^. 

l<j<k<N j=l 

Assume a G {-1,1 + 1}, b G {-l /2,{l + l)/2}, h G {-h/2,{h + l)/2}, b 2 G 
{-/ 2 /2, (Z 2 + l)/2}, and 6 3 e {-Z 3 /2, (Z 3 + l)/2}. 
As k — > while a, 6 05 &i> ^2, b 3 are fixed, 

(4.4) (-e(x) o Y x o G^)- 1 + C )/n 2 -> tf, 



where if is the Hamiltonian of the -BCjy Inozemtsev model given in (|2.1| ) and Co is 
a constant. Hence, we recover the Hamiltonian of the BCn Inozemtsev model from 
a operator of BCn Ruijsenaars model via a limit k — > 0. 

Let us make a correspondence between the invariant spaces of theta functions on 
the BCn Ruijsenaars model and the space related to the quasi-exact solvability on 
the BCn Inozemtsev model. 

Proposition 4.2. Let a = — fi/n, bo = —{vi + Pi)/2k, b\ = —[y 2 + v 2 )/2k, b 2 = 
— (i/s+Ps)/2k, 63 = — (z/q+^o)/2k. LetTh^} BN ^ ( U--Q ) be the theta-type invariant space 
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of BCn Ruijsenaars model, Wj7 m ( 3.11 ) be the invariant space of BCn Inozemtsev 
xlel, and Q(x) 

+ 63) e z> . 

Then the map 



•> " V - 7 K u- 1 / " - "J -J " 

model, andQ(x) be the function defined in ( \j-4 )- Assume k = —((N — l)a + b + bi + 
b 2 + 63) e Z> . 



(4.5) 0: ThZ (BN) - W* ym 

f(x!,...x N ) 1 ► 6(x)/(xi, . . .Xjv) 

zs an isomorphism of vector spaces. 

Proof. Let us consider the correspondence between the space ThZ (BN) and the space 
V fc sym , where Vj! ym was defined in proposition pTT| . 

Let /(xi, . . . ,2^) G Th^ k (BN) andgfa, ...,x N )= f(x t , . . . ,x N )Q(x)$(p(xi), p{x N 
where the function . . . , z^) was defined in ( |2.2|) . From the condition fc £ Z, the 
function ^(xi, . . . , £jv) does not have branchs on C^. It is seen that the function 
g(xi, . . . ,xn) is doubly periodic, W(B N )-mva.na.iat, and may have poles only along 
Xj = (J = 1, . . . , N) up to periods with degree at most k. 

Hence there exists g(zi, . . . , zn) G V fc sym such that g(p(xi), . . . , p(xjsr)) = g(x\, . . . , xn) 



by a similar argument in the proof of proposition pT2. 



By composing with the canonical map from V fc sym to W^ ym } we obtain (j)(Th^ BN ^) C 
W^ ym . It is obvious that the map is injective, and the dimension of Th^ Bjv ^ is 
equal to that of Wc fm . Therefore the map <fi is bijective. □ 



In proposition |4~2| , we have established that the theta-type invariant space Th^ k 
of the BCn Ruijsenaars model corresponds to the space W^ ym which is related to the 
quasi-exact solvability of the BCn Inozemtsev model. 



5. Degenerate Inozemtsev model 
5.1. Trigonometric BCn Inozemtsev model. 



In ||16|| , Sasaki and Takasaki considered degenerate BCn Inozemtsev models and 
their quasi-exact solvability. They also considered for the case of type An- 

In this section, we consider the degeneration of BCn Inozemtsev model and show 
that the finite-dimensional invariant spaces for the elliptic BCn Inozemtsev model 
tend to the spaces introduced by Sasaki and Takasaki by the degeneration. 

The Hamiltonian of the trigonometric (or degenerate) BCn Inozemtsev model is 
given as follows: 

N 



j=l i l<j<k<N ^ 



7T 2 7T 2 



+ 



sin 2 ir(xj — Xk) sin 2 n(xj + Xk) 



N 



+£ 



* M ° + ^ + ^i'Lt ^ + Sccanr, + See, to, ) , 



Sin 7TXj COS^ TTXj 



where I, lo, h, ci, C2 are coupling constants. 

This model is known to be quantum integrable. In [13|], Oshima gave the explicit 



expression of commuting operators of conserved quantities. 



12 



KOUICHI TAKEMURA 



Set 
(5.2) 



$ D (x) 



- N 

a 



N 



exp j — — y~^cos 2nXj J J^[(si 



Sill 



\/ +i, 



COS TXX, 



\h+l 



3=1 



j [ (sin 7r (x^ — Xj 2 ) sin tt (xj 1 + Xj 2 ) ) 
i<ii<i2<iv 



l+i 



2m 2 



Let be the vector space spanned by functions &d(x) (sin 7rxi) 2mi (sin 7rx 2 * 
such that rrii G {0, 1, ... , L} for all i, and H / | D '*' sym be the set of ©^-invariant ele- 



(sm ttxn 



,2m/> 



ments in W 



(D) 



The following proposition is essentially shown in ]T6| . 

Proposition 5.1. (c./. [16) The Hamiltonian ( \5.J\ ) preserves the space W^ ,sym , 
ifLe Z> 0; c 2 = and ~c x = 2~aix 2 (2L + l + h + 3 + 2{N - 1)(/ + 1)). 

Proof. We set Wo = log<&o(x). Then we have 



« S((S 



+ 



<9 2 iy 

dx 2 



E 

l<j<k<N 



27T 2 /(/ 



+ 



2tt 2 /(/ + 1) 



sin 7r(x 7 - — Xfc) sm 7r(a; 7 - + 



N 

+ E 

i=i 



7r 2 / (/ + l) TT^xgH-l) 

. 9 ~r n 



Sin TTXj 



COS^ 7TX 



C3 cos 2iiXj 



vr 2 a 2 



COs47TXo 

2 3 



+ C 



where C is a constant term and c 3 = 2a7r 2 (l + li + 3 + 2( iV — 1) (I + 1) ) . By comparing 
with the Hamiltonian in (7.1)] and its corresponding 'exactly solvable sector' fT6| , 
(7.13)], we obtain the proposition. □ 

5.2. Degeneration. 

In this section, we consider the degeneration (the trigonometric limit) r — > \J — loo 
and see the correspondences of Hamiltonians and their invariant spaces between the 
nondegenerate model and the degenerate one. 

Let I, Iq, h, I2, h be the coupling constants of elliptic Inozemtsev model (see (|2.1|) ). 
We adopt the following limits of coupling constants as r — ► \J — loo: 

• I, l , l\\ fixed; 

• h — ^V 1 + b an d ^3 = — fP -1 + b\ where p = exp(7r-\/— It). Here we note that 
p —>■ as r — > \J — loo. Then the Hamiltonian if of the elliptic Inozemtsev model 
(see ( |2.1| )) tends to the Hamiltonian of the trigonometric Inozemtsev model (see 
(p.ip). More precisely, 

2 N 

(5.4) 



2 JV q2 

^+y(^2 + l) + / 3 (/3 + l))--E^2 



+ E 

l<j'<fc<JV 



2tt 2 /(/ 



+ 



Sin 7Tl 



2vr 2 /(/ + 1) 
sin 2 it (xj + Xk) 



N 

E 



vr 2 / (/o + I] 



Sin TTXj 



n 2 h(h + l) 



COS^ TTXn 



2vr 2 fi(2&+ 1) COS27TXo 



7T 2 S 2 



COS 47TXj I + C; 



as p — > 0, where Ci is a constant. 
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Let us observe how the invariant space varies as p — > 0. Since b = - 2 ^ 2 - 1 the 
Hamiltonian H (see ( |2.1|) ) preserves the space W^ m (see ( p.lll) ) for a — I + l,&i = 
^±i, 6 2 = =f, and 6 3 = ^ if L = -(N - I) (I + 1) - + 6 e Z> . 

We consider the limit p -»• 0. Note that if L = -(iV - 1)(Z + 1) - k±|i±g + g e ^> 
then the Hamiltonian if preserves the space W^ m whenever p varies. 

Let $(z) = Ui< j<k <N( z 3 ~ Zk) l+1 Uj=ifa ~ e i)^~( z j ~ e 2 )^(^ - e 3 )^ be 
the function defined in Q for a = Z + l,6i = = ^,& 3 = Then 

$(p(xi), . . . , p(xjv)) — > C3^_d(^) as r — > y/^-loo, where C 3 is a constant and 



N 



(5.5) * D (z) = n(sin7r^)- 2(iV - 1)(Z+1) - (Zl+1)+2b (cos7rx 




[ [ (sin ^(xjj — Xj 2 ) sia7r(xj 1 + Xj 2 )) /+1 exp 

l<Ji<j2<AT 

Set t(x) = sin ^ — — Let W L be the vector space spanned by functions ^ d{%) 
t( Xl ) mi t(x 2 ) m2 . . • t(x N ) mN such that rrn e {0, 1, ... , L} for all z, and H>^ D) ' sym be the 
set of ©Ar-invariant elements in 

As p — > 0, the vector space H^ ym tends to the space Wf° ^' sym ) and the operator 
which appears on the right-hand side of (|5.4j ) preserves the space W^ ,sym if L = 
—(N - 1)(Z + 1) - + 6 G Z> . If Z + 1 = 2b - 2(N - 1)(Z + 1) - (Zi + 1) then 

it is seen that W^ ,sym = W^ ,sym . Therefore, we recover proposition |5.1| by the 
trigonometric limit. 

In summary, by the trigonometric limit we have shown that some finite-dimensional 
invariant space of the Hamiltonian H of elliptic model (see (|2.1)) tends to the the 
invariant space of Sasaki and Takasaki which is related to the quasi-exact solvablity. 

In [I3|, Oshima described the limit procedure of the commuting operators of con- 
served quantities. By applying Oshima's result and Proposition |3]4] in this article, it 
follows that the commuting operators of the trigonometric BCn Inozemtsev model 
also preserve the space W^ ,sym . 

Hence we established that the commuting operators of the trigonometric BCn 
Inozemtsev model also preserve the space related to the quasi-exact solvability. 



6. Concluding remarks 

In this paper, quasi-exact solvability for the BCn Inozemtsev model is proved not 
only for the Hamiltonian but also for commuting operators of conserved quantities. It 
is seen that the theta-type invariant spaces for the BCn Ruijsenaars model correspond 
to the spaces which are related to the quasi-exact solvability for the BCn Inozemtsev 
model, and the degeneration of the BCn Inozemtsev model (especially for its quasi- 
exact solvability) is clarified. 

In papers |^, EO], several models which are related to the Inozemtsev model are 
studied. It would be interesting to link their results with ours. In [ [TD] , I8|, the method 



of perturbation for the elliptic Calogero-Moser-Sutherland models from the trigono- 
metric models is introduced. For the Hamiltonian of the BCn Inozemtsev model, 
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holomorphy of perturbation from the trigonometric model can be established. Rela- 
tionship between the perturbation and the complete integrability should be clarified. 
More precisely, holomorphy of perturbation for commuting operators of conserved 
quantities should be shown, although it is not successful as of now. 



Appendix A. 

This appendix presents the definitions of and formulae for elliptic functions. 
Let uj\ and u 3 be complex numbers such that the value uj^/uji is an element of the 
upper half plane. The Weierstrass p-function is defined as follows: 

(A.l) p(x) = p(x\%j u 2u)s) = 

J_ / 1 1 \ 

x 2 , s ^-^ „ \ (x — 2muJi — 2nu 3 ) 2 {2muJi + 2noos) 2 J ' 

(m,n)eZxZ\{(0,0)} VV 6J V 7 

Setting cu 2 = — ct> 3 and 

(A.2) ei = p(ui) (i = 1,2,3). 

yields the relations 

(A.3) ei + e 2 + e 3 = 0, p(x + 2uj) = p(x); (j = 1, 2, 3), 

1 ( 1 + 1 + 1 



p'{x) 2 2\p(x)-e 1 p(x)-e 2 p{x) - e 3 J ' 

^ + y) + p(,-,) = gg±^-2p(,)-2p( y ), 

p(x + ^) = e t + (e -~7y et ~ et//) , (< = 1,2,3), 

where i', i" G {1, 2, 3} with i' < i", i ^ i', and i ^ i" . 

Let wi = 1/2 and t — uo^ / uo\ . The Jacobi theta functions are defined by 



(A.4) e x {x) = 2 Y^{-l) n - l e™^ n -^> 2 sin(2n - l)irx, 

n=l 

oo 

e 2 ( x ) = 2J2 e T7r ^ n -^ 2 cos(2n - l)nx, 

n=l 

oo 

3 (x) = 1 + 2 e TUv/3In2 cos 2mix, 

n=l 

oo 

6 (x) = 6i{x) = 1 + 2^(-l) n e W ^ Tn2 cos2n7ra; 



n=l 
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Then the following relations are satisfied: 

d 2 

(A.5) p(x) = -—log (x) + const, 

d 2 

p(x + u)i) = -— log0 i+1 (x) + const (i = 1, 2, 3), 

= -Oi{-x), 8i(x) = 6i{-x) (i = 0, 2, 3), 
O^x + 1) = -fl^x) (i = 1, 2), Oiix + 1) = 0,(x) (i = 0, 3), 
^(x + r) = - e -v^T(2x+r)^ (;r) ( - = 0> 1)? 

^(x + r) = e-^^+^x) = 2, 3), 
0i(2x)0 2 (O)0 3 (O)0 o (O) = 20 1 (x)9 2 (x)9 3 (x)9 o (x), 

d / 1 (o) = nd 2 (o)e 3 (o)e (o). 



Let p = exp(7Tv^— It). The expansions of the functions p(x), p(x + |), p(x + |) 
and p(x + i^ 1 ) in p are given as follows: 

2 2 00 2w 

7T 7T Z 2 ^ np 



(A.6) p(x) = , 2 - - - Stt 2 J2 T-^( cos2 ^ " !)■ 
sin ra 3 1 — 7r n 

v ' n=l 
/ 1 \ ^2 2 °° 2n 

p ( * + o = — ^ - t - 87r2 £ r^^- 1 )" cos2mnr - x )> 

V 2/ cos (7rx) 3 1 — p 2n 

\ / \ > n=l 

/ t\ 7T 2 9 ^ cos 2,irnx — p ra 

H* + 2) = -¥- 8 ' r I> ■ 

n=l ^ 

/ l + r\ 7T 2 2 ^ „(-l) n cos27rnx-p n 
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